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We report on the first observation of Di(2420) — > Dtv + it~ decays (where the contribution from 
the dominant known D\ — > D*tx decay mode is excluded) in the B — > Din decays. The observation 
is based on 15.2 x 10 7 BB events collected with the Belle detector at the KEKB collider. We also 
set 90% confidence level upper limits for the — > D'*V + 7t~ and D\ — > D*it + tt~ decays in the 
B -» £>2(2460)7r and B -» Di(2420)tt decays, respectively. 

PACS numbers: 13.25.Hw, 14.40.Lb 



The D± — > D*n decay is currently known to be the pri- 
mary decay mode of the Di(2420) meson 0. However, 
the transitions D± — > Dtttt via other intermediate quasi- 
two-body resonance states or via non-resonant decays are 
possible and may contribute to the D\ total width. Mea- 
surements of their branching ratios and analysis of the 
decay dynamics are particularly relevant to a study of 
production rates of various D** excitations in B decays. 

The ratio of the branching fractions R = B(B~ — > 
D*°tt-)/B(B- -> £>?7r-) is calculated in HQET and the 
factorization approximation in Refs. 0,0]. InRef. 0, R 
is found to depend on the values of the subleading Isgur- 
Wise functions (ti j2 ) describing h.QCD/m c corrections; 
thus measurement of R can be used to estimate the sub- 
leading functions. In Ref. 0, some of the subleading 
terms are estimated and the ratio is determined to be 

(l+*? a )/(l+*? 1 )' 



R w 0.35 



where J? 1 ^ 2 ^ are non- 



factorizable corrections that are expected to be small. 

The first observation of D** (denotes P-wave D excita- 
tions) production in B decays was reported by CLEO . 
From their studies and the measurement of the ratio 
B(D*° -» D+tt-)/B{D*° -» D*+ir~) the R value 
was determined to be R — 1.8 ± 0.8, where it was as- 
sumed that decays of Di and D 2 mesons are saturated 
by the two-body Dir, D*ir modes. Recently the branch- 
ing fractions for the decays B — > D**n — > D^mr have 
been measured with better accuracy 0, resulting in 
R = 0.77 ± 0.15. The existence of Dn >2 ) decay channels 
other than D\ — » D(*'tt can also affect the R value, either 
decreasing or increasing the currently observed 2.8ct dif- 
ference between the prediction and experimental results. 
In this Letter we report the first observation of the 

D°ir-ir+ decays. The D x 
and 



D+ -> D+7T-TT+ and D\ 
mesons were reconstructed from the B 



o 



DiTT decays, respectively. The results are based 



on a sample of 15.2 x 10 7 BB pairs produced at the KEKB 
asymmetric energy e + e _ collider The inclusion of 
charge conjugate states is implicit throughout this report. 

The Belle detector has been described elsewhere 0. 
Charged tracks are selected with a set of requirements 
based on the average number of hits in the central drift 
chamber (CDC) and on the distance of the closest ap- 
proach to the interaction point. Track momentum trans- 
verse to the beam axis of at least 0.05 GeV/c is re- 
quired for all tracks in order to reduce the combinatorial 
background. For charged particle identification (PID), 
the combined information from specific ionization in the 
CDC (dE/dx), time-of-flight scintillation counters and 
aerogel Cerenkov counters is used. Charged kaons arc 
selected with PID criteria that have an efficiency of 88%, 
a pion misidcntification probability of 8%, and negligi- 
ble contamination from protons. All charged tracks with 
PID responses consistent with a pion hypothesis that are 
not positively identified as electrons are considered as 
pion candidates. Photon candidates are selected from 
calorimeter showers not associated with charged tracks. 
An energy deposition of at least 30 MeV and a photon- 
like shape are required for each candidate. Pairs of pho- 
tons with an invariant mass within 12 MeV/c 2 (~ 2.5<r) 
of the 7T° nominal mass are considered as ir° candi- 
dates. 

We reconstruct D°(D + ) mesons in the K~ir + 
(K~tt + tt + ) decay channel and require the invariant 
mass to be within 15 MeV/c 2 (~ 3cr) of the D°(D+) 
mass. Then, D*°(D* + ) mesons are reconstructed in the 
D°ir°(D°ir + ) decay mode. The calculated mass differ- 
ence between D*°(D* + ) and D° candidates is required 
to be within 2 (1.5) MeV/c 2 (~ 2.5a) of the expected 
value 0. For D* -> D°tt decays the D° -> i^^Tr+Tr+vr" 
mode is also included (the same D* parameters were used 
as above). 
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TABLE I: Number of events, efficiencies and branching fraction products of B — » D**tt,D** — > D'*'7r + 7r decays. 



Mode 


iV sl9 


e (10- 2 ) 


B (10~ 4 ) 


Significance 


B~ -> Dln~ 


, D? -► L>°7r-7r+ 


151 ± 24 


14.1 


(1.85 ±0.29 ±0.35^ 


-0.0 \ 
-0.46 1 


8.7a 


B° D+7T- 


, L>^ -> D + 7T-7T+ 


124 ± 20 


9.9 


(0.89 ±0.15 ±0.17j 


-0.0 \ 
-0.26 1 


10a 


fi- -> D?7T 


, D? -> ZTV+tt" 


< 1.2 


2.2 


< 0.06 






5° -» D+7T- 


,D+ -> D* + TT + TV- 


< 12.0 


3.4 


< 0.33 










< 4.4 


2.2 


< 0.22 






B° -> DJ+tt 


~,0; + -> D* + 7r + 7r- 


< 9.0 


3.4 


< 0.24 







We combine candidates with 7r~7r~7r + to 

form B mesons. Candidate events are identi- 
fied by their center-of-mass (CM) energy difference, 
AE = Q2i Ei) — fibcam, and the beam constrained mass, 
M bc = V^boam _ (12iPi) 2 > where E hcam is the beam 
energy and Pi and Ei are the momenta and energies of 
the decay products of the B meson in the CM frame. 
We define the signal region as 5.273 GeV/c 2 < M bc < 
5.285 GeV/c 2 and |A£| < 25 MeV. The sidebands are 
defined as 5.273 GeV/c 2 < M bc < 5.285 GeV/c 2 and 
25 MeV< \AE\ < 50 MeV. If there is more than one B 
candidate in an event, the one with Z)W mass closest to 
the nominal value and the best 7r~7r-7r + vertex is cho- 
sen. We use Monte Carlo (MC) simulation to model the 
detector response and determine the acceptance 10]. 

Variables that characterize the event topology calcu- 
lated in the CM frame are used to suppress background 
from the two-jet-like e + e~ — ► qq continuum process. We 
require | cos#thr| < 0.80, where 9 t hi is the angle between 
the thrust axis of the B candidate and that of the rest 
of the event; this eliminates 77% of the continuum back- 
ground while retaining 78% of the signal events. We also 
define a Fisher discriminant, T, which is based on the 
production angle of the B candidate, the angle of the 
thrust axis with respect to the beam axis, and nine pa- 
rameters that characterize the momentum flow in the 
event 11]. We impose a requirement on T that rejects 
67% of the remaining continuum background and retains 
83% of the signal. 

To suppress the large contribution from the dominant 
D\ — ► D*ir — ► Dtttt decay mode we apply a requirement 
on the invariant mass of the relevant Dtt combination 
\{m Dv -m D )-{m^ G -m^ DG )\ > 6 MeV/c 2 l] (10a). 

The AE and M D (.)„ distributions for the selected 
B — > D±Tt, Di — > D^ww candidates are shown in 
Fig. To plot the AE distributions, we require Mb c 
to lie in the signal region with an additional require- 
ment \M D (.)„ n - M Dl \ < 25 MeV/c 2 , where M Dl is 
the Di world average mass value; for the Md(.) w dis- 
tributions we select events from the AE signal region. 
(Although there are two Dtt + tt~ combinations, they are 
kinematically separated in the D\ mass region.) Clear 
signals are observed for B~ — > D^n^ , D® — > D Q ir~ir + 
and B° — > D\tt~ , D\ — > D + tt~'k + decays. For branch- 



ing fraction calculations we use signal yields determined 
from the fit to Momt distributions as it allows us to di- 
rectly estimate a possible contribution from B — > Z?27r, 
Di — > Dtttt decay. The signal shape distribution is pa- 
rameterized by a convolution of a resolution Gaussian 
(a = 2.5 MeV/c 2 ) with a signal Breit-Wigner func- 
tion; the background is represented by a linear func- 
tion. The D\ mass and width determined from the fit are 
M D o = 2426 ± 3 ± 1 MeV/c 2 (statistical and systematic 
error, respectively), T D o = 24±7±8 MeV/c 2 for and 
M D + = 2421 ± 2 ± 1 MeV/c 2 , T D + = 21 ± 5 ± 8 MeV/c 2 

for Dt; these are consistent with the world average val- 
ues |l|. The signal yields are given in Table Q (in the 
cited branching ratios the first and second errors are sta- 
tistical and systematic; where a third error given is due 
to model uncertainty). For the B — ► D\-k — > £>*7r~7r~7r + 
decay channels, we do not observe statistically signifi- 
cantly signals and thus determine 90% CL upper lim- 
its |l2j for their branching fractions. In the fit to the 
MjjtnTT distribution, we fix the D\ mass and width at 
their world average values. The statistical significance of 
signals quoted in Table [I] is defined as yj— 2 In (Lo / L max ) , 
where L max and Lq denote the maximum likelihood with 
the nominal signal yield and with the signal yield fixed 
at zero, respectively. 

To account for contamination from other possible D\ 
production mechanisms (such as e + e~ — > cc continuum 
production or semileptonic B — > D\lv decays), we fit the 
Mb-ktt distribution for events in the AE sidebands. In 
this fit, we fix the D\ mass and width at their world 
average values. The fits give —6 ± 8 events for the D® 
and 10 ± 11 events for the D^. 

The B -> D|(2460)tt,LI2 -> Dtt+tt" decay may 
also contribute to the B — > Z?7r~7r~7r + final state. 
To analyse a possible effect, we made a simultaneous 
fit of the M(£)?7T + 7r _ ) and M{D% vr+vr") distributions, 
where we assume isospin invariance and require the ratio 
iV(D|)/iV(£)i) to be the same for both charge combina- 
tions. The fit finds the ratio N{D* 2 )/N{D 1 ) = 0.33±0.14 
and signal yields of N(D\) = 120 ± 17, N(D+) = 
107 ± 16. Thus, we set a 90% CL upper limit for D*: 
B(B -> D* 2 tt-) x B[D* 2 Dir+ir-) < 0.55B{B -> 
Di7T _ ) x B{Dx -> Dtt+tt') 13]. The number of events 
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for the D\ mechanism obtained in the fit with the -pp- 
ratio fixed to 0.46 gives the model uncertainty for the D\ 
yield: j^i^- Another asymmetric uncertainty, coming 
from the other possible D\ sources, is tw% for D° and 
^22% for D±- These two uncertainties are combined in 
the final resutls as "model" uncertainty. 

The signal yields extracted from the AE distributions 
are used only for a consistency check of the results. The 
AE signal shape is parameterized by a Gaussian with 
parameters determined from signal MC. The AE back- 
ground shape is described by a linear function. We re- 
strict the fit to the range -0.1 GeV < AE < 0.2 GeV to 
avoid contributions from other B decays, where an addi- 
tional pion is not reconstructed. Signal yields obtained 
from the fits to AE distributions are 106±12 for D°tt + 7t~ 
and 96 ± 13 for D + tt + tt~ , while the corresponding recon- 
struction efficiencies are 10.8% and 7.6%, respectively. 
Thus the obtained reconstructed event numbers by the 
two methods are consistent within the statistical uncer- 
tainty. 

In order to determine the total D\ — > Dtttt width, anal- 
ysis of final states with neutral pions is required. With 
only the D\ — > Dtt + tt~ branching fraction measurement, 
the analysis of the decay dynamics could also be use- 
ful to determine the total D\ — > Dtttt width. As the 
limited statistics do not allow us to perform the full am- 
plitude analysis, we consider the one-dimensional projec- 
tions of several variables: M D „, M^,-, cos Q^^tt^,, ), 
cos9(7r i3 7r^„), and cos8(tt^-D) (where all angles are 
calculated in the D** rest frame). Although these vari- 
ables are not independent, they highlight each model's 
features. For instance, the helicity angle distributions 
differentiate between the D\ — > D(tttt) and D\ — > 
(Dtt)tt models. We select events from the B sig- 
nal region with the additional requirement \M (Dtttt) — 
Md x I < 25 MeV/c 2 . Decays through the following quasi- 
two-body intermediate states are considered: D\ — > 
Dp -> Dtt+tt-, Di -> Z^(2308)tt -► Dtttt and D x -> 
D/ o (600) -> Dtt+tt- (we set M fo = 0.8 GeV/c 2 and 
r /o = 0.8 GeV/c 2 ; the 1^(2308) parameters are taken 
from Ref. 0). We use the simplest non-trivial Lorentz- 
invariant expressions for the corresponding matrix ele- 
ments in MC simulation 01 . We fit the experimental 
data with different models. For each variable we plot two 
distributions: one from the signal region and the other 
from the AE sideband. We perform a simultaneous fit 
to these distributions, assuming a Poisson-like profile in 
each bin whose mean is the sum of the background and 
signal (for a given model) in the signal region or the back- 
ground only in the sideband. The obtained differences of 
likelihood values for all variables are listed in Table ITT1 
Figure |21 shows the M„+ v - and cos<d(TTgD) distributions 
along with expectations based on different D\ — > Dtt + tt~ 
decay models 0. Although the D\ — > DqTt decay mech- 
anism describes the data best, some contribution from 



TABLE II: Comparison of goodness-of-fit tests for the con- 
sidered D\ — > Drnr decay models. 



Distribution 


{-2\nL Dp /L D> f/ 2 


(~2\n LofJLo^) 1 / 2 




D° , D+ 


D° , D+ 


M D7I 


2.5 , 2.7 


1.9 , 2.9 




3.4 , 1.6 


5.2 , 4.7 


COS 0(7Tg7rJ»») 


1.6 , 2.5 


-2.0° , 3.4 


COS 0(7Tg7T^*» ) 


2.5 , 3.0 


4.0 , 2.9 


cos ©(n^D) 


2.0 , 0.5 


3.2 , 4.0 



"In this case Lj>j < L D »^ 



other mechanisms cannot be excluded completely. 

It is interesting to examine the dependence of the R 
value on the decay mechanism. The expression for R can 
be written as s x x B{B~ -> Dfn-^Df -> D*+tt-)/ 
(si x B(B~ -> Dlir~,Dl -> D*+tt-) + s 2 x B(B~ -» 
DItt^jD® — > Z?°7r + 7r - )) where Si is a scale factor that 
recovers the full width from the single decay channel. 
(It includes branching fractions of other possible de- 
cays of both D** and a meson from its decay prod- 
ucts.) Following the procedure used in Ref. and 
fixing si at 3/2, we can calculate s 2 factors for differ- 
ent models: S2(D**° — > Dq + tt~) = 9/4 (disregarding 
possible interference effects in D\ — > D + tt~tt° decays), 
s 2 (D**° -> D°p) = 3, s 2 {D**° -> £> /oL= 3 / 2 - Using 
the branching fractions measured in Ref. [jj and here, the 
central value for R depends on the decay model: 0.50 for 
Dp, 0.60 for Df and 0.54 for D^tt. 

The following sources of systematic errors are con- 
sidered: tracking efficiency (8% overall, integrated over 
particle momenta), kaon identification efficiency (2% 
overall), tt° reconstruction efficiency (8%), D branch- 
ing fraction uncertainties (2%-7%), MC statistics (2%), 
model uncertainty in MC efficiency (10%), uncertainty 
caused by variation of cuts (5%), background shape un- 
certainty (10%). The uncertainty in the tracking effi- 
ciency is estimated using partially reconstructed D* + — > 
D°[Kc ! tt + tt~]7t + decays. The kaon identification uncer- 
tainty is determined from D* + — » D q [K~tt + ]tt + decays. 
The 7T° reconstruction uncertainty is obtained using D° 
decays to K~tt + and K~tt + tt°. To determine the system- 
atic uncertainty in the signal yield extraction, we use dif- 
ferent parameterizations for the background events. The 
overall systematic uncertainty is 19% for B — » Dtttttt 
and 21% for B — > D*tttttt. We assume equal production 
rates for B + B~ and B°B° pairs and do not include the 
corresponding uncertainty in the total systematic error. 

The B~ — ► D°tt + tt~tt~ final state also includes the 
D* + tt~tt~ intermediate state with D* + — > D°tt + . We re- 
verse the D* veto requirement to select D* + tt + tt~ events 
and measure the branching ratio B(B~ — > D* + tt~tt~) — 
(1.27 ± 0.07) x 10~ 4 (based on a sample of 85 x 10 6 BB 
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events), that agrees well with the value of B(B — > 
D*+ir-ir-) = (1.25 ± 0.07) x 10~ 4 measured earlier 0. 

In summary, we report the first observation of 
Z?i(2420) — > £>7r + 7r~ decays (with the dominant Di — » 
Z)*7r contribution excluded). The measured branch- 
ing ratios with the corresponding statistical significances 
and systematic uncertainties are presented in Table [i] 
We find the upper limit for the possible D\ contri- 
bution to these results: B(B — > D^n^) x B{D* 2 — > 
Dit+t:-) < 0.556(5 -> £>i?r-) x B(Z>! -> Dtt+tt-). 
No statistically significant signal has been observed for 
the D** — > L>*7r+7r~ decays. The corresponding 90% 
CL upper limits are listed in Table [I] Analysis of the 
D\ — ► Z?7r + 7r~ dynamics shows that the decay model 
Di — > Dg7r gives the best description of the data. The 
R = B(B~ -> Dfn-)/B(B- -> DJtt") value calculated 
assuming Di — > -Dp 71 " dominates is 0.54 ± 0.18; this is 
~ 2a lower than the previously published one. 

We thank the KEKB group for the excellent operation 
of the accelerator, the KEK Cryogenics group for the ef- 
ficient operation of the solenoid, and the KEK computer 
group and the Nil for valuable computing and Super- 
SINET network support. We acknowledge support from 
MEXT and JSPS (Japan); ARC and DEST (Australia); 
NSFC (contract No. 10175071, China); DST (India); the 
BK21 program of MOEHRD and the CHEP SRC pro- 
gram of KOSEF (Korea); KBN (contract No. 2P03B 
01324, Poland); MIST (Russia); MESS (Slovenia); NSC 
and MOE (Taiwan); and DOE (USA). 



* on leave from Nova Gorica Polytechnic, Nova Gorica 
[1] S. Eidelman et al., Phys. Lett. B 592, 1 (2004). 
[2] A.K. Leibovich, Z. Ligeti, I.W. Stewart, M.B. Wise, 

Phys. Rev. D 57, 308 (1997). 
[3] M. Neubert, Phys. Lett. B 418, 173 (1998). 
[4] M.S. Alam et al. (CLEO Collaboration), Phys. Rev. D 
50, 43 (1994). 

[5] P. Avery et al. (CLEO Collaboration), Phys. Lett. B 331, 
236 (1994) [Erratum-ibid. B 342, 453 (1995)]. 

[6] H. Albrecht et al. (ARGUS Collaboration), Phys. Lett. 
B 232, 398 (1989). 

[7] K. Abe et al. (Belle Collaboration), Phys. Rev. D 69 
112002 (2004). 

[8] S. Kurokawa and E. Kikutani, Nucl. Instr. and Meth. A 
499, 1 (2003). 

[9] A. Abashian et al. (Belle Collaboration), Nucl. Instr. and 

Meth. A 479, 117 (2002). 
[10] R. Brun et al., GEANT 3.21, CERN DD/EE/84-1, 1984. 
[11] D.M. Asner et al. (CLEO Collaboration), Phys. Rev. D 

53, 1039 (1996). 
[12] G.J. Feldman and R.D. Cousins. Phys. Rev. D 57, 3873 

(1998). 

[13] As most of the systematic uncertainties cancel in the ratio 
of the branching fractions, the systematic uncertainty is 



negligible. 

[14] We use the following expressions for matrix elements for 




AE (GeV) M(D*ji*7t") (GeV/c 2 ) 




AE (GeV) M(D**7i*7t") (GeV/c 2 ) 



FIG. 1: AE (left) and Md™ (right) distributions for the 
Di — > D°-k~-k + (first row), Df — > D + tt~-k + (second row), 
D° -» D* ( \~it + (third row) D+ -> D* + 7r"7r + (fourth row). 
Open histograms represent the data from the signal area, 
hatched histograms show the Md-wtv (where applicable) and 
AE sidebands, respectively, the curves are the fit results - for 
the signal area and sidebands. 



the Di -» Dtvtt decay: F^^B , ir+)F^(Do, 7r~_ )■ 

P{p[\P{pzY, F^(B°,7r+)F^(f,D-y 
P(D^) v P(D^r, F^(B°,n+)F^(p°,n+ - 
P(D 1 ) V P(D 1 ) tl , where the following notation is used: 
F liV {A,B) = P>XP V B - P V A P£, P(B - C) = P(B) - P(C) 
and P stands for a 4-momentum. 
[15] In case of the D fo model the key cos Q(ngD) distribution 
is practically independent of the fo mass and width. 
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FIG. 2: M n+7T - (left) and cose(7r~D) (right) distributions 
for the Di — > D°-k~-k + and Df — » D + 7r _ 7r + , respectively. 
Points with error bars represent the experimental data, solid 
line - D*°iv, dashed - Dp, chain - Dfo models with the ex- 
pected background added. The hatched histogram corre- 
sponds to expected background (from AE sidebands). 



